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Abstract 



The hydrodynamic limit for the Boltzmann equation is studied in the case when 
the limit system, that is, the system of Euler equations contains contact discon- 
tinuities. When suitable initial data is chosen to avoid the initial layer, we prove 
^ I that there exists a unique solution to the Boltzmann equation globally in time for 

any given Knudsen number. And this family of solutions converge to the local 
00 . Maxwellian defined by the contact discontinuity of the Euler equations uniformly 

away from the discontinuity as the Knudsen number e tends to zero. The proof is 
'sj" ■ obtained by an appropriately chosen scaling and the energy method through the 

. micro-macro decomposition. 

o 

, J . 1 Introduction 

' Consider the Botlzmann equation with slab symmetry 

ft + U. = ^Q{fJ), (/,a;,t,0eRxRxR+xR3, (1.1) 

where ^ = (^i, ^2, ^3) ^ R^, f{x, t,C,) is the density distribution function of the particles at 
time t and space x with velocity ^, and e > is the Knudsen number which is proportional 
to the mean free path. 

The equation (11.11) was established by Boltzmann [1] in 1872 to describe the motion 
of rarefied gases and it is a fundamental equation in statistics physics. For monatomic 
gas, the rotational invariance of the particles leads to the following bilinear form for the 
collision operator 



1 



2 



F-M. Huang, Y. Wang, T. Yang 



where ^',^1 are the velocities after an elastic collision of two particles with velocities 
before the collision. Here, 6 is the angle between the relative velocity ^ — ^* and the unit 
vector in = {f2 G : — ■ > 0}. The conservation of momentum and energy 
gives the following relation between the velocities before and after collision: 

e = e-[(e-e*)-fi] fi, 
e = e* + [(e-e*)-fi] fi. 

In this paper, we consider the Boltzmann equation for the two basic models, that is, 
the hard sphere model and the hard potential including Maxwellian molecules under the 
assumption of angular cut-off. That is, we assume that the collision kernel B{\^ — C,*\, 0) 
takes one of the following two forms. 



and 



B{\i-^ 



B{\i-i. 



n — D _ 
n — 1 i 



Here, n is the index in the inverse power potential which is proportional to r^^" with r 
being the distance between two particles. 

Formally, when the Knudsen number e tends to zero, the limit of the Boltzmann 
equation fll.ip is the classical system of Euler equations 



Pt + ipui)x = 0, 
{pui)t + {pul +p)x = 0, 
{pui)t + {puiUi)cc = 0, i 
\u 



where 



[piE + 

p{x,t) = 
pui{x,t) 
P{E + 



2,3, 



'1.21 



+ [pui{E + — ) +pui 



R3 



1,2,3, 



u\ 



-)ix,t) 



;i.3) 



Here, p is the density, u = {ui,U2, us) is the macroscopic velocity, E is the internal energy 
and the pressure p = Rp9 with R being the gas constant. The temperature 9 is related 
to the internal energy hy E = lR9, and <^i{0{i = 0,1,2,3,4) are the collision invariants 
given by 

' MO = 1, 

MO = for ^ = 1,2,3, 

MO = iiep, 



that satisfy 



/ MOQih,9)d^ = 0, for 2 = 0,1,2,3,4. 
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How to justify the above limit, that is, the Euler equation (11. 2p from Boltzmann 
equation (11.11) when Knudsen number tends to zero is an open problem going way back 
to the time of Maxwell. For this, Hilbert introduced the famous Hilbert expansion to 
show formally that the first order approximation of the Boltzmann equation gives the 
Euler equations. On the other hand, it is important to verify this limit process rigorously 
in mathematics. For the case when the Euler equation has smooth solutions, the zero 
Knudsen number limit of the Boltzmann equation has been studied even in the case with 
an initial layer, cf. Asona-Ukai[r|, Caflishp], Lachowicz|21j and Nishida[29] etc. However, 
as is well-known, solutions of the Euler equation (11.21) in general develop singularities, such 
as shock waves and contact discontinuities. Therefore, how to verify the hydrodynamic 
limit from Boltzmann equation to the Euler equations with basic wave patterns is an 
natural problem. In this direction, Yu [35] showed that when the solution of the Euler 
equation (II. 2p contains non-interacting shocks, there exists a sequence of solutions to 
the Boltzmann equation that converge to the local Maxwellian defined by the solution 
of the Euler equation (11.21) uniformly away from the shock. In this work, the inner and 
outer expansions developed by Goodman-Xin [12] for conservation laws and the Hilbert 
expansion were crucially used. 

The main purpose of this paper is to study the hydrodynamic limit of the Boltzmann 
equation when the corresponding Euler equation contains contact discontinuities. More 
precisely, given a solution of the Euler equation (II. 2p with contact discontinuities, we 
will show that there exists a family of solutions to the Boltzmann equation that converge 
to a local Maxwellian defined by the Euler solution uniformly away from the contact 
discontinuity as e — >■ 0. Moreover, a uniform convergence rate in e is also given. The proof 
is obtained by a scaling transformation of the independent variables and the perturbation 
together with the energy method introduced by Liu- Yang- Yu|24] . 

For later use, we now briefly introduce the micro-macro decomposition around the 
local Maxwellian defined by the solution to the Boltzmann equation, cf. [2l]. For a 
solution f{x,t,^) of the Boltzmann equation (II. ip . we decompose it into 

f{x,t,0=M{x,t,0 + G{x,t,0, 

where the local Maxwellian M(x, t,^) = M[p„5i](^) represents the macroscopic (fluid) 
component of the solution, which is naturally defined by the five conserved quantities, i.e., 
the mass density p{x,t), the momentum pu{x,t), and the total energy p{E + ^\u\'^){x,t) 
in (ll.3p . through 

^'"^^ ' v/(27ri?^(x,t))3 ^ ^ 

And G(x,t,^) being the difference between the solution and the above local Maxwellian 
represents the microscopic (non-fiuid) component. 

For convenience, we denote the inner product of h and g in L|(R^) with respect to a 
given Maxwellian M by: 

7r3 M 
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If M is the local Maxwellian M defined in (11.41) . with respect to the corresponding inner 
product, the macroscopic space is spanned by the following five pairwise orthogonal base 

Xo({) = -^M, 
f- 

{XhXj) = Sij, hi = 0, 1,2,3,4. 

In the following, if M is the local Maxwellian M, we just use the simplified notation (-, ■) 
to denote the inner product (■,-)m- We can now define the macroscopic projection Pq 
and microscopic projection Pi as follows 



1.5) 



^oh = ^{h,Xj)Xj, 

j=0 

Pih = h-Poh. 
The projections Pq and Pi are orthogonal and satisfy 

PoPo = Po, PiPi = Pi, PoPi = PlPo = 0. 
We remark that a function h{^) is called microscopic or non- fluid if 

I hiOfiiOd^ = 0, 2 = 0,1,2,3,4, 

where is the collision invariants. 

Under the above micro-macro decomposition, the solution f[x,t,C,) of the Boltzmann 
equation (11. ip satisfies 

Po/ = M, Pi/ = G, 
and the Boltzmann equation (11.11) becomes 



(M + G)t + ei(M + G), = -[2g(M, G) + Q(G, G)]. 



1.6) 



If we multiply the equation (II. 6p by the colhsion invariants (pi{^){i = 0,1,2,3,4) and 
integrate the resulting equations with respect to C, over R'^, then we can get the following 
fluid-type system for the fluid components: 



pt + [pui]^ = 0, 
{pui)t + (puj+p), 

{pui)t + {puiUi)x = 



^i^iGr^d^, i = 2,3, 



1.7) 



[p(E + ^-^)], + [pmiE + ^-^) + pu,l = - j -m'G.dC 
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Note that the above fluid-type system is not closed and we need one more equation for 
the non-fluid component G which can be obtained by applying the projection operator 
Pi to the equation ( 11.6^ : 

Gt + Pi(6M,) + Pi(6G,) = i [LmG + Q(G, G)] . (1.8) 

Here Lm is the linearized collision operator of /) with respect to the local Maxwellian 
M: 

^mh = 2Q(M, h) = Q(M, h) + Q{h, M). 
And the null space M of Lm is spanned by the macroscopic variables: 

J = 0,1,2,3,4. 

Furthermore, there exists a positive constant ao{p,u,6) > such that for any function 

MOeAA^, cf. 

< h, liMh >< -o-Q < jy{\^\)h, h >, 

where I'd^l) is the collision frequency. For the hard sphere and the hard potential with 
angular cut-off, the collision frequency has the following property 

o<z^o<Kl^l)<c(l + lel)^ 

for some positive constants z/q, c and < k < 1. 

Consequently, the linearized collision operator Lm is a dissipative operator on L^(R^), 
and its inverse L^^ exists and is a bounded operator in L^(R^). 

It follows from ([L8]) that 

G = eL^[P^{^^M,)] + e, (1.9) 

with 

e = L^[e{Gt + Pi(6G.)) - Q(G, G)]. (1.10) 
Plugging the equation fll.9p into (11.71) gives 

pt + (pui)x = 0, 

f 



{pUi)t + (pMlMi)x' = £ipiO)uix)x - j ^iCiQxd^, « = 2, 3, 

I M P I M P As 

[p{E + '-^)]t + [pui{E + ^) + pui]x = e{\{e)ex)x + —{p{e)uiui, 
+eY^{p{d)uiUix)x - I ^^il^Pe^c^^, 



1.111 



where the viscosity coefficient p{9) > and the heat conductivity coefficient A(^^) > are 

2 
3 



smooth functions of the temperature 9, and we normalize the gas constant i? to be | so 
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that E = 9 and p = \pd. The exphcit formula of p{9) and \{9) can be found for example 
in [36], we omit it here for brevity. 

Since the problem considered in this paper is one dimensional in the space variable 
a; G R, in the macroscopic level, it is more convenient to rewrite the equation f ll.ip and 
the system (11.21) in the Lagrangian coordinates as in the study of conservation laws. That 
is, set the coordinate transformation: 



X 



p{y,t)dy, 



t t. 



We will still denote the Lagrangian coordinates by (x, t) for simplicity of notation. Then 
(11. ip and (II. 2p in the Lagrangian coordinates become, respectively, 

V V e 



and 



( Vt- Ui^ 

Ult + Px 
Uit 

'u 



0, 
0, 



0, z = 2,3, 



+ 



+ {pu 



0. 



Also, ffLTD- ffTTTj) take the form 



Vt - Mix 
Ult + Px 



e^Gxd^, 



Gt 



with 



Uit = - ^i^iGxd^, i = 2,3, 
ImP r 1 

+ + (m). = - J ^m'Gxd^, 

^G, + -Pi(eiM,) + iPi(eiG,) = -(LmG + g(G, G)), 
V V V e 

G = £LM'(ipi(eiM,)) + ei, 

Bi = L^[e{Gt - ^G, + ^Pi(eiG,)) - Q{G, G)]. 



and 



0, 

ie^fi{9) 

y 



Vt - 



V 



,p{9) 



-Ui 



U.<S>i.d^, ^ = 2,3, 



+ 



\u\ 



2 -)t + iPUi)x = £{ 

( UiUi^)^ ■ 

V 

1=2 



m 



4e.fi{9) 



Ox)x + -77- 

V 3 



-UiUix 



;i.i2) 



;i.i3) 



;i.i4) 



;i.i5) 

;i.i6) 
;i.i7) 



:i.i8i 
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In the following sections, we will apply some scaling and energy method to these 
equations. 



2 Main result 

We will state the main result in this section. For this, we firstly recall the construction 
of the contact wave {v,u,6){x,t) for the Boltzmann equation in ^18j. Consider the Euler 
system ( I1.13P with a Riemann initial data 

(-•^"-°'-{|::;o:::i; :tl <-> 

where v±,6± are positive constant. It is well-known (cf. [30j) that the Riemann problem 
(ll.lSp . (12. ip admits a contact discontinuity solution 

(V^,[/,e)(x,t) = | (2.2) 

provided that 

p. := = p+ := . 2.3 

Motivated by (12. 2p and (12. Sp . we expect that for the contact wave {v,u,9){x,t), 

m 

V 

Then the leading order of the energy equation (ll.lSP /i is 



p = — ~ p+, <^ 1. 



0,+p^u^, = e{^^),. (2.4) 

V 

By using the mass equation (ll.lSp i and w ~ we obtain the following nonlinear diffusion 
equation 

et = 6iaie)9.U aie) = ^^^. (2.5) 

From [2], |9], we know that the nonlinear diffusion equation (12.50 admits a unique self- 
similar solution Q(r]), rj = . ^ with the following boundary conditions 

Ve(l+*) 

0(-oo, t) = 0(+oo, t) = e+. 

Let 6 = — 6^\. Q{x,t) has the property 

ejx,t) = , ^ e as x->±oo, (2.6) 
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with some positive constant c depending only on 6±. 
Now the contact wave {v,u,9){x,t) can be defined by 



-6, Ml 



2ea{Q) 



3p+ 3p 
Note that {v,u,9){x,t) satisfies the following system 
Vt - Uix = 0, 

Ult+Px = —{^—Ulxjx + Rlx, 

-3 V 

Uit — cl — z — 



0, (i = 2,3), e = e- 



\u\ 



+ 



\u\ 



+ {pu 



2,3, 



l)x 



V 3 V 



where 



Ri 



+£( > , -^UiUixjx + n.2x 



-^a(e)0i + (p - J9+) - ^^u,x = 0(l)fe(l + t)-^e-^, 

i?2 = -(A(e)e. - x{9%) + {p- p+Yh - ^^u^u,, 

V „ 3v 



with some positive constant c > depending only on 6±. 
From (I2.6p . we have 



0{l)6e 2^(1+*), 
0(l)(5e"2^^, 



Therefore, 

' \{v,u,e){x,t) - {v_,0,e_)\ = 0(l)5e"WTiy 
\{v,u,9){x,t) - K,0,e+)| = 0{l)6e-m 
We are now ready to state the main result as follows. 



if X < 0, 
if X > 0. 



if X < 0, 
if X > 0. 



(2.7) 



(2.8) 



(2.9) 
(2.10) 



:2.11) 



(2.12) 



Theorem 2.1 Given a contact discontinuity solution {V,U,Q){x,t) of the Euler sys- 
tem U.13\) . there exists small positive constants 6o, Eq and a global Maxwellian M^, = 
^[p*,u,,e:t], such that if 5 < 5q, e < Sq, then the Boltzmann equation ( li.i)) admits a unique 
global solution f'^{x,t,^) satisfying 



R3 



\f'{x,t,0-M[y^tjM^,t,0\' 



■d^ < Ci5oe^ + 02606 



(2.13) 
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with some positive constants Ci{i = 1,2,3) independent of e. 
Consequently, we have 

sup \\r{x,t,i) - Myy^tj,e]{^,t,mm^) <Ch5oe-^, V/i > 0, (2.14) 

\x\>h «^vT^^ 

where the norm \\ ■ ||j-2( ^i ^ is ||--7=||^2(j^3). 

Remark. Theorem 2.1 shows that, away from the contact discontinuity located at x = 
0, for any Knudsen number e, there exists a unique global solution /'^(x,t,^) of the 
Boltzmann equation (11.11) which tends to Mjy f; 0](a:, t, ^) as two global Maxwellians with 
a jump at X = when e ^ Moreover, a uniform convergence rate in the norm 
^^^Ktm:) holds. 

3 Reformulated system 

In this section, we will reformulate the system and introduce a scaling for the independent 
variable and the perturbation. Firstly, we define the scaled independent variables 

y = e ^x, r = e ^t. (3.1) 

Correspondingly, set the scaled perturbation as 

v{x,t) = v{x,t) +e^(f){y,T), 
u{x, t) = u{x, t) + e^ip{y, r), 

e{x,t) = e{x,t)+eH{y,r), 

{e + ^)ix,t) = {e + Lj-){x,t)+e-^uj{y,T), 
G(x,t,0 = ^^G(y,r,0, 

ei(x,t,o = £^ei(y,r,o. 

We remark that the above scaling transformation plays an important role in the following 
proof. 

Under this scaling, the hydrodynamic limit problem is now transferred into a scaled 
time-asymptotic stability of the viscous contact wave to the Boltzmann equation. In fact, 
this scaling is suitable for the contact wave because of its parabolic structure. Notice that 
the hydrodynamic limit proved by this method is globally in time unlike the case with 
shock profile proved in [35] which is locally in time. However, we do not know whether 
there exists some appropriate scaling for the shock profile so that this method can be 
applied. 

With the above scaling, the proof of Theorem 2.1 will be given by energy method as 
[T8] for the scaled perturbation (0, C)(|/, r) and G{y,T,^). 

From the construction of the contact wave {v,u,d), the relation between the viscous 
contact wave {v, u, 9) to the Boltzmann equation and the inviscid contact discontinuity 
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(V^, U, G) is given by fl2.12p . Thus, in order to prove Theorem 2.1, it is sufficient to consider 
the convergence of the solution /(y, t, ^) of the Boltzmann equation to the Maxwelhan 
M[^^^ g](?/, r, ^) defined by the contact wave {v,u,9) as the Knudsen number e tends to 
zero. 

For this, as in \18\, we introduce the following anti-derivative of the perturbation: 



(3.3) 



Obviously, 

From fll.lSp and (12.81) . we have the following system for ($, \E', W) 



- ^ly = 0, 



3' V 



-U. 



a^.Gide, z = 2,3. 



Wr + £ 2 [pui — pui) 
3 



.x{e)^ \{e)^^ 4,/i(e) ii{e) 



'y _ "yl 

V V 



3' V 



UiUiy - e 2 



i=2 



(3.4) 



where the error terms Ri (z = 1, 2) are given in (12. 9p and (12.101) . 
Introduce a new variable 

W = W -ui-^i. 



(3.5) 



It follows that 



C = Wy-Y, with Y = -e--\^!y\'' -Uly^l. 



(3.6) 



By using the new variable W and linearizing the system (13. 4p . we have 



' - ^^y = 0, 

f _ 3f 3 f 



1) 



1 

£2- 



wy 



Wr+p+^iy = e-^^Wyy- j^m'^id^ + ui jeiQid^ + Q^, 



(3.7) 
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where 



3 V 3v 

Q4 = ^ - ^^^j, + —^-^uiy'^iy -e 2i?2 - Mir^i + e"uiRi 
V V 3 V ^ " 

1=2 



and 



(3.9) 



We now derive the equation for the scaled non-fluid component G(y, r, ^). From (11.150 . 
we have 

- ^G, + e~-2-P,{^,My) + -Pi(6G,) = £-5LmG + Q(G, G). (3.10) 

V V V 

Thus, we obtain 

G = ^Lj^^[Pi(eiM,)]+ei, (3.11) 

and 

e,{y, T, = e'^L^lGr - ^G, + ^Pi(eiG,) - Q{G, G)]. (3.12) 



Go(?/,r,0 = ,^L^{P^U^-^^Oy + ^■Uy)M]}, (3.13) 



Let 

^'"'"-"'^ = 2.^"----^ 29 
and 

G^{y, T, = Giy, r, - Goiy, r, 0- (3.14) 

Then Gi{y,T,^) satisfies 

G.-.-iLMG,= -J^P,|6(^C. + «-«Ml ^^^^^^ 

+^G, - iPiKiG,) + Q(G, G) - Go,. 

Notice that in (13.141) and (13.151) . Go is substracted from G because WOyW^ ~ (l + e^r)"^/^ 
is not integrable with respect to r. 

Finally, from (I1.12p and the scaling (13.11) . we have 



fr--fy + -fy = e-'^Q{f,f). 



ly < Jy-^ -^\J,J J- (3.16) 
V V 
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In the following, we will derive the energy estimate on the scaled Boltzmann equation 
(I3.16p . Indeed, to prove Theorem 2.1, it is sufficient to prove the following theorem. 

Theorem 3.1. There exist small positive constants 6i, ei and a global Maxwellian M* = 
^lv,,u,,e*] such that if the initial data and the wave strength 6 satisfy 

^6(r)|.=o + 5<5i, (3.17) 

and the Knudsen number e satisfies e < £i, then the problem Ii3.16\) admits a unique global 
solution f'^{y,T,^) satisfying 

sup||r(i/,r,0 -M[^^^^e-](?/,r,0llL2( 1 ) < C6ie^. (3.18) 

y 

Here Eq{t) will be defined in ( 15. ip satisfying 

e,{t)^ m,^,w)r+\\{<i>,^x)r+e\\{K^yXy)r+ J J^/^dy 

^ f f\d'''G\\^^ ^ f f\d"f\\,, * (3-19) 

From now on, 9", d°'' denote the derivatives with respect to y or r, and || ■ |p represents 
II ■ III2 for simplicity of notations. 

Remark: In particular, if we choose the initial value of the Boltzmann equation fl3.16l) 
as 

r(y,0,O = M[,,^,e-](|/,0,O = M[,(^,o),^M),%,o)](0, (3.20) 

then 

E6(r)U = 0(l)[||(e"„n,)f + e||(,;,„^,„^,,)fl =0(1)5. (3.21) 

In fact, the initial data f{y, 0, ^) can be chosen such that the initial perturbation Eq{t)\t-=q 
is suitably small and of order 0(1) with respect to e. This is the reason why we use the 
scaled variables y,T in (13. ip . otherwise, the initial perturbation EQ{T)\r=o is not uniform 
with respect to e. 



4 A priori estimate 

We will focus on the reformulated system (13.71) and (I3.15p . Since the local existence of the 
solution to (13. 7p and (13.150 is now standard, cf. [31] or [36], to prove the global existence, 
we only need to close the following a priori estimate by the continuity argument 

N{T)= sup l\\{<^,^,w)\\i^ + \\{<p,^,or+em,^y,Cy)r+ f f ^« 



0<T<T 

-e 



\a'\=l'^ |a|=2 



(4.1) 
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where 7 is a small positive constant depending on the initial data and the strength of the 
contact wave, and M^, is a global Maxwellian chosen later. 

We now briefly explain the a priori assumption ||($, \E', W^)||icx) < 7^ in (14.11) . Roughly 
speaking, based on the observation in \l8l that the energy estimate involving || (<l>, \E', W) ||^2 
may grow at a rate (l + ear)^, the decay of \\{^x, ^x, W^a;)|li2 in the order of (l + £2r)~2 is 
needed to compensate this growth. This yields a uniform boundedness of ||($, VI^)||ioo, 
which is essential to close the a priori estimate. 

Note that the a priori assumption (14. ip also gives 

^^II(0,V^,C)IIL<^7^ (4.2) 
_ _ 1. _ 1 

^^11 / g<«IU. < Csi (/ / grf?*,)' . (/ / ^-^didyj < + (4,3) 

and for \a'\ = 1, 

.ill /^..u. . CI (/ p^^^il /M.,.,)* , a,..r 

(4.4) 

From (ICTl and (ESD, we have 



= -e 2 {—Uiy)y - / iliiGydi, « = 2, 3 



Cr + £ ' [VUly - pUly) = -e '^[^—Oy)y —[^—UiUiy)y-e ^ R2y 

1) O V 



i,\u\^ 



V 

Thus 

ell ((/.., ^.,C.)f<C(5 + 7)'. (4.6) 

Hence, we have 

II {v^. 9,) f < Cell (0., 7^., C.) IP + C^ll {vr, 9r) f < C{5 + 7)'. (4.7) 
In addition, (14.11) also implies that 

II {Vy, Uy, Qy) f < CS\\ (0„ ^y, (y) ^ + C \\ (Vy, Uy , Q y) f < C {5 + 7)^. (4.8) 

Since 

e\\d-(p,pu,p{E+''^)\\\'<CEj p.^didy<C^\ (4.9) 
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(ra-iSD give 

e\\d^{v,u,e)f <Ce\\d^[p,pu,p{E + 



\u\ 



+CeY: I |9"'(p,P«,p(i^ + ^))rrfl/ (4.10) 

Thus, for |a| = 2, we have 

^, f < Ce{\\d^{v, n, e)f + u, ) < C(5 + 7)'. (4.11) 
Finally, from the fact that / = M + e^G, we can obtain for |a| = 2, 



< 



Ce I l\^d^dy + Ce\\d"{v,u,e)r + Ce J \d-' iv,u,e)\'dy (4-12) 

\a'\=l 



l«'l = 

<C(5 + 7)'- 

Before proving the a priori estimate fl4.1l) . we list some basic lemmas based on the cele- 
brated H-theorem for later use. The first lemma is from [13] . 

Lemma 4.1. There exists a positive constant C such that 



/ 



M u M J m J m J M 

where M can be any Maxwellian so that the above integrals are well defined. 

Based on Lemma 4.1, the following three lemmas are proved in [25]. The proofs are 
straightforward by using Cauchy inequality. 

Lemma 4.2. If ^/2 < < then there exist two positive constants a = cr(f , u, 6; f^,, u^,, 9^ 
and ?7o = riQ{v,u,9;v^,u^,9^) such that if |f — v^\ + \u — -u*! + \9 — 9^\ < rjQ, we have for 

Lemma 4.3. Under the assumptions in Lemma 4.2, we have for each h{^) & A/"^, 

-^\-L^h\di<a J ^^^e, 

m^-M^d,<.-^[^-^d,. 
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Lemma 4.4. Under the conditions in Lemma 4.2, for any positive constants k and A, it 
holds that 

where the constant Ck^x depends on k and A. 
4.1 Lower order estimate 

Now we will derive the lower order estimates of ($, \E', W). By multiplying fl3.7l) i by 
( I3.7P 9 by v'^i, f l3.7P ci by ^'j, {\3.7\i a by respectively and adding all the resulting 

equations, we have 

1=2 ^+ 1=2 ^+ 

i=2 ^+ 

2W 

+vQi^!^ + 2^ Qi^!, + — Q4 + iVi + (■ ■ Oy, 

i=2 ^+ 

(4.13) 

where 

iVi = -^^i j il^idi-Y^^^ j U^Qldi+'^{u^ j il^idi- j \ii\i?QidO- (4-14) 

From now on, (■ ■ ■)y denotes the term in the conservative form so that it vanishes after 
integration with respect to y over R. Let 

We estimate the right hand side of (14.131) term by term as follows. Firstly, 

jhr^ldy<C6e-2{l + eh)-'Ei, (4.16) 

/ ^He)]y^i'i^iydy<PK, + CpSeHl + eh)-'Eu (4.17) 
where /3 is a small positive constant to be chosen later. 



(4.15) 
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Now we estimate f vQi^idy by 

+ j \e-^^vRi^i\dy + j \^Y^i\dy:=^Ii. 



(4.18) 



(4.19) 



(4.20) 



Note that 

h<Ce"^ j \{%,Ou^y^My + Ce j \{%X)^iy'^i\dy 

h <C j {e^^y\^ + e'^lWyl'' + eW^ + \u\'')\^<i\dy 

< C6e^l + ehy^Ei + C-f{e^\\^yf + Ki) + C6e^l + eh) 

h<C5£^l+eh)-^Ei + C6e^l + eh)-^, (4.21) 

and 

/4<Cy" \{e^\^y\^ + uiy<ili)<ili\dy <C56^1 + eh)-^Ei + C-fKi. (4.22) 
Substituting fHl8D - (l4:22D into KWl yields 

J vQi^idy < Cfc^(l+£^r)-^Ei+C(5+7)(£^||<l>J|2+iri)+Cfci(l+£^r)-^+C7£i||7/'iJ|2 

(4.23) 

Similarly, we can estimate 

f f 2W 

J Qi^idy {i = 2, 3) and J -^Q^dy. 

Now we estimate / Nidy. We only need to estimate Ti =: — J -u^i J ^fOid^dy because 
other terms in j Nidy can be estimated similarly. Let M,, be a global Maxwellian with 
the state (f*, u*, 9^,) satisfying < 6^ < 6 and \v — \ + \u — u^ \ + \6 — < tiq so that 
Lemma 4.2 holds. By the definition of Gi, cf. (13.121) . we have 

Ti =-£^ j j ei1^MiGr)d^dy + e'^ j j eii^MiGy)d^dy 



— £2 

4 



e,L^[P,i^,GyMdy + e'^ / v^^ / e?Lj^i[Q(G, G)]^^^^ (4.24) 



i=l 



For the integral T^, we have 

TI = -e^ j vm, j eiL^{Gir)didy - e'^ j j ei^m{Gor)didy 
=:Tl^+Tl\ 



(4.25) 
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Note that the hnearized operator satisfies, for any h E A/""*", 



(4 26) 

L^h)y = L^{hy)-2L^{Q{L^h,My)}. 



Then we have 



-2e-2 j v^, j eii^M{Qi^MGi,Mr)}d^dy. 



(4.27) 

J e?L^^(GOrfer <cj "-^\G,\'d^. (4.28) 



The Holder inequahty and Lemma 4.3 yield 



Moreover, from Lemmas 4.1-4.3, we have 



I e?Lj^HQ(LMGi,M.)K<c(y" '^\L^{Q(L^G„M^)}\'d^ 
Combining (OTjl - fH:^ gives 



:5 j j eil^^Gid^dy)r + C5e-^{l + e"^r)-^E^ 



On the other hand, by fl3.13p . we have 

Tr = -e"- j v^i j ei'i^MiGor)d^dy 

<Ce5 / \^iJ^\{\{9yr,Uyr)\ + \iey,Uy)\\iVr,Ur,er)\)dy 

< C6e^{l + e"^T)-'E^ + C6e{l + ehyl + C6el \\ (0,, C.) f , 
which, together with fl4.30p . imply 

Tl < -{e^ j v^i j ^lL^Gid^dy)r + C5e^ (1 + ehy^Ei + Cfe(l + eh) 



(4.30) 



(4.31) 



(4.32) 
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The estimation on Tf, i = 2,4 is relatively easy by using the Cauchy inequality and 
Lemmas 4.1-4.3. In fact, direct computation gives 



Tt<Ce-^ I I '^\Gy\'d^dy + Ce'2 j mluldy 

<C5e^{l + e^T)-^E^ + C^eK^ + Ce"^ [ [ ^^\Gy\^d^dy. 




On the other hand. 



Ti < Cie'^ I ( / ^|Lm {Q(G, m'dO'^dy 

Klel),^,2, 




Then, it follows from (3.46), (3.55) and Lemmas 3.1-3.4 that 

_ 4 



-2^^ /^/ ^l'i^M{Qi'i^M[Pii^iG)],My)}d^dy 



(4.33) 



<C^e-^ J J -^\G\'d^dy (AM) 
<C7£5 [ f '^Ml\G,\''d^dy + C6e^l + eh)-K 



The estimation on Tf is similar to the one for T/. Firstly, notice that 

4 

Pi(eiG,) = [Pi(eiG)], + ^ < eiG, X, > Pi(x,g. (4.35) 



< C6eHl + e^ry^Ei + C(7 + (3)K^ + Cfe^ (i + 52^)- 2 + C^ei || (0^^, (j^)||2 

(4.36) 



+Ce^^ I I ^\G,\'d^dy. 




By IKm . gSZD-dlSID and fOB]) . we have 



+CPe'^\\^,4' + C{j + P){K, + e^\%r) + Ce'^ j j ^Ml\G,\^didy 



-Ce^ I I ^\Gy\'d^dy + Ci5 + ^)sl ||a"'(0,^,C)f 

\a'\=l 

(4.37) 
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The estimates on the other terms of J Nidy are similar and we omit the details for 
brevity. Therefore, collecting the above inequalities gives 




$, VI/, W)l.^{G^)didy)r + < Cife5 (1 + e"^r)-^E, 

|2 , , .A,i|i^ ii2 , A f f '^(l^l)ir^ |2, 



1 , i_ 

" 2 



|a'|=l 



/ / '^\Gy\'d^dy + C^i5 + ^)el ^, Of + + e^r) 

(4.38) 

where we have used the smallness of 5, (3 and 7. Here W) is a linear function 

of ($, ly) which is a polynomial of ^. 

Note that the dissipation term Ki does not contain the term e5||<|)^||2. To complete 
the lower order inequality, we have to estimate From (8.8)2, we have 

l£lM^$ - + ^-±% = ^Wy -Q,+ f eiQidt (4.39) 
3 V V 3v J 

Multiplying (14.391) by e^^y yields 



(4.40) 



' 3v 

where we have used 

$,^lr = ($,^l).-($.^l), + ^?,. 

Integrating fl4.4Up with respect to y gives 

< CKi + Ce^2 j Qldy + Ce^ j \ j ^fB^d^l'^dy. 
By (13.81) and the Cauchy inequality, one has 



(4.41) 



2||„/, ||2 

lyW ■ 



j Qldy < C6e^^il+6'^T)-^Ei+Cje{Ki+e'^\\%f)+C6e^l+eh)-'2+C{6+-f)e 

(4.42) 

On the other hand. Lemmas 4.1-4.3 imply 



I ' / ^CelJ2 J J '^Id'-'Gl'd^dy 



l"'l=i (4.43) 
+C-fe j j ^^|Giprferft/ + Cfei(l+£5r)"i. 
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Thus combining (032])- (033]) yields 

< C2&4(1 + Eirj-'E, + CjA'i + C2&4(1 + eiryi + C(S + 7)£^l|i;'i!,f (4.44) 

\a'\=l 

The microscopic component Gi can be estimated through the equation (13.151) . Mul- 
tiplying p.l5p by ^ gives 

29 

+^G, - -Pi(eiG,) + Q(G, G) - Go.} 



Gi 



Integrating fl4.45p with respect to ^ and ?/ and using the Cauchy inequality and Lemma 
4.1-4.4 yield that 




|a'|=l 



On the other hand, from the fluid- type system (13. 7p . we can get an estimate for 
£3 11(^^,1^^)11 2 as follows. 



3 , 



£5 II (^,, Wr) f < C^e-- (1 + e^T)-'Ei + C^K^ + C^e^^ W^W" + C4£5 1| (^^, 
+C45£i(l + £5r)^i + C4((5 + 7)£^ /" 



|a'|=l 



M, ' ^' (4.47) 



We can now complete the lower order estimate. Since A(^, PF) is a linear function 
of the vector ($, 14^) which is a polynomial of ^, we get 

iL^'Girferfyl <]E, + Ce [ [ ^d^dy. 




We choose large constants Ci > 1, C2 > 1, C3 > 1 and small constant P such that 




C-iEi + Cie^^ I j AL^Gid^dy + C2 j e-^^^^^ '^^ 



E2 = CiEi + CiE-^ I I AL^Gid^dy + C2 I e^^<!>l-s-^%^^dy 



3v 



+C3 / / ^didy (4.48) 
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and 

Hence, multiplying (OS]) by Ci, flCTI) by C2, (OB]) by C3, (I07|) by Ci{6 + 7)^1 and 
adding all these inequalities imply 



^ " ~7 , (4.50) 

\a'\=l 

where 

i^2 = ^i^i + C72|e^|±$Jrfy + .^||(M/.,iy.)f + ^^3.-^ I l^\G,\'d^dy. (4.51) 
4.2 Higher order estimate 

In this subsection, we shall estimate the derivatives of Applying dy to the 

system (13.41) gives 



(pT - 'iply = 0, 

6 V V J 



J ^l^^Qlyd^, ^ = 2,3, 



\(f)\ \(0) 

Cr + e^^ipUiy - pUly) = ( 0y — ^ y) y + Q5 

3 



1=1 

where 

3 



i=2 



(4.53) 



Multiplying (14.521) 9 by ipi, (I4.52P '^ by ipi {i = 2, 3) respectively and adding them together 
yield 

3 — 



1=1 

3 ,n\ /n\ 3 



+ J^( ^ij; —Uiy)'4)iy = -e'^Riytpi - X^^i / ^l^i^lyd^ + 

i=2 V V J 
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Since p — p = £2 + R0{1 — i), we obtain 



V V V 6 V " ^-^ V ^ 

1=1 1=2 



3 



A V V i=i J 



(4.54) 



Set 

<l(s) = s - 1 - Ins. 

Then 

{Re^-)}r = -e^ROi- - - p^i-)vr + (4.55) 

V V V V V 

where 

^(s) = s"^ - 1 + Ins. 

It is easy to check that $(1) = <I>'(1) = ^(1) = ^'(1) = and $(s) is strictly convex 
around s = 1. Substituting fl4.55p into fl4.54p yields 

3 < 3 



^-^ 2 V V 6 V ^-^ V V 

i=l 1=2 

— 3 

+e-^vpM-) - ^(^ - ^)uiy^l)iy-e''^Riy^i -y2^^ f ^i^i^iyd^+i 

V 6 V V ^ J 



Note that 



(4.56) 

{ml)}r = eHl-^-)Cr-i>{l)er. (4.57) 



and 



(4.58) 
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Substituting fH37D and (gSHD into fH36|) gives 



3 



,i=l / r 1=2 _ 

"U f U o V V 



(4.59) 

where 



Let 



3 3 

I ^l^^Qlydi+^^{Y.U, I i,i,Q,ydi-]- I m^Qlydi). (4.60) 

i=l ^ i=l ^ ^ 



3 

^ (4 61) 

i=2 

Integrating (I4.59P with respect to y yields 

E-,, + ^K^<C6ekl + eh)-'\\{%,^y,Wy)f + C6e-2{l + e"^T)-l+ j N^dy. (4.62) 

Here, we only consider the term — / V^i / C,fQiydC,dy because other terms in J N2dy can 
be estimated similarly. By fl4.43p . one has 

~ I I ^1®^^*^^^^ 
i^iy / il^id^dy 



<\K, + C6eHl + e'^r)-UCe^ J2 [ [ "-^{d^' G\^d^dy ^^'^^^ 



\a'\=l 



+C{5 + ^) J J '^IG.l'd^dy. 
Combining IK62\i ) and IK63\i yields 

\a'\=l * * 
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We need to estimate £2 which is not contained in K^. Following the same way as 
in estimating £2||$^||2 in the previous subsection, we firstly rewrite the equation (14.521) 9 
as 

\^^^-^^yT - i'lr - ^~^{P - P)y 



(J O U J 



3 

by using the equation of conservation of the mass (14. 521) i . 
Since 

-{p-p)y = e^-(py - e^—Cy + ( -)vy - ( —)0y, 

" V V V V V V 

and 

then by multiplying (I4.65P by e^cpy, we get 

i'-^ec^l - eHyi.,). + e\^l = ei^Wy + ehly + .^f 

-(^ - l)My + i§ - §Wy<Py - l<^)yi'iy<Py (4-66) 



'Uly]y<Py + Rly4>y + j ^l^lyd^^y 



3 V V 

Integrating (I4.66P with respect to y and using the Cauchy inequality yield 

< CjKs + CjSe^l + e^r)-i || {<i>y, ^y, Wy) f + Cj6el (1 + ehyl + C^-fe 



|a|=2- ' 

Here we have used 



(4.67) 



e'^ [ I / e^elyd^\'dy < Cel [ [ ^\d^G\'d^dy 



+C^(* + ^>E / l'^\^'G\'di''y^C(S + -,)ei [ ['MGA^ii,y (4.68) 




\a'\=V 

-C6e\e^\\(j)yf + K^) + Cfet(l + e^r)"!. 



To estimate £2 1|(0^, -i/;^, ^^)||2^ we need to use the equation (14.51) . By multiplying (l4.5P i 
by £20^^ (I4.5P 9 by e^ipir, (14.50 -? by £2^.^ (^ = 2,3) and (l4.5P /i by £2^^ respectively, and 
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adding them together, after integrating with respect to we have 




Klel),^ ,2^.^., (4.69) 



Choose large constants C4, C5 > 1 such that 



and 



Let 



E, = C,E, + C,j {^e<j>l - e^y^,)dy + j j ^d^dy, (4.70) 

i^4 = ^i^3 + f /e^|0Sc^y + 5^||(0.,V^.,C.)f + /^|Gir«, (4.71) 

then from (ICTD . (14:671) and flOTj) . we have 

+^9.^ E / / ^|5"'GP« + C9.t 1 1 '-^\d-G?d^dy. 

|a'|=l |c»l=2 

(4.72) 

Next we derive the estimate on the higher order derivatives. By multiplying (14.521) 9 
by —eipiyy, (14. 521) c! by —eipiyy {i = 2,3), (I4.52P /L by s(yy, and adding them together, we 
obtain 

i=l i=2 

4 , , A sMie), , , 3,A(e), , , 

- (— — JyV'iy^iyy - 2^ {——)y^iy^iyy ' {——)yCyCyy 
o (J '2 

-^^[^^ - ^)^l,]y^l,. - ^[(^7^ - -^)dy]yCyy + sHp - P)y^iyy (4-73) 
o (J (J (JO 

+62 Riyljj^yy + e^{puiy - pUiy)Cyy - eQ^Cyy + ^'^i'iyy / 

1=1 

3 

where Q5 is defined in (14.531) . 
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Integrating (14.731) with respect to y yields 



i=l " i=2 

|2 I 7>- \ I r^r^i/i I iTr TJ/ M|2 



3 



<C(e'||.#.,||- + A3) + C&i(l+£'T)-'||(*„*„»y||" + C&i(l+£'r)-3 

3 I 



101=2- 

Now we get the estimation of By applying dy to (14 .5 1) 9. we get 



i'lyr + e ^^{p-p)yy = -^{^^uiy)yy-e Riyy - j ^iGyyd^. (4.75) 



Note that 



- = + -Cj/y --yP- P)^yy ' ^ ' -Pm --f{p-P)y- -f<Py (4-76) 

Multiplying fITTSD by —£'2(j)yy and using (I4.76P imply 



-{j e-^lPly(t)yydy)r + ^ £ ^ _ ^^^rf?/ < ^ 1 1 (^^^^^ , (^^^ ) 1 1 ^ + ^6 ^ ( 1 + £ 2 T ) 2 

+C5e^(l + 5^r)-l($„ VI/,, W^,)f + C(5 + 7)(^^ll0.f + K,) (4.77) 



I"l=2 



To estimate e'i\\{(t)yrjipyTXyT)\\'^ and £2 1| (0^^, -^Z;^^, ^^^) ||2^ we use the system (14. 5 p again. 

2 3 

By applying dy to (14.50 . and multiplying the four equations of (14.51) by e^(t)yTi ^'^i'lyr, 

3 3 

^'^i'iyr {i = 2,3), e2(y^ respectively, then adding them together and integrating with 
respect to y give 



£^ II (0?/r, V'j/r , Cyr) f < C^^^ || ^j/y, Cyj/) f + C6e^l + sh) i 
+Cfe^(l + £^r)-i(<l>„ VI/,, iy,)f + C{6 + 7)(£^||0,||' + Ks) 

+Ce^ J J "^IGyl^d^dy + Cel j j "-^W^Mdy. 

\a\=2 



(4.78) 



3 

£2 



Similarly, we have 

rT,ijTrXrT)f < || (0,^ , (,^) || 2 + ( 1 + t) " t 

+C6e^l + e'^r)-^\\{%,my,Wy)f + C{5 + ^)e'^ ^ ||a"'(0, ^, OH' 

|a'|=l (4.79) 

\a\=2 
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By choosing Cq and C-j to be large enough, we have 

3 „/,2 /•2 



|q!|=2 |q!'|=1 

1 _L ch 'T-\~'^ W f (h \T/ U 

3/ 



+C(5 + 7M E l|5"'(0,^,C)f + Cfe^(l + ^^r)-i||($„vl/^,iy,)f 

|a'|=l^ 

To close the a priori estimate, we also need to estimate the derivatives on the non-fluid 
component G, i.e., d^G, {\a\ = 1,2). Applying dy on (I3.15p . we have 



e~-^LMGy + 2E'-^Q{My, G) + 2QiGy, G). 



Since 



we have 



|{-Pi(6M,)},| < Civ'y + u'y + e'y + \eyy\ + |n,,|)|i?(0|M, 



yy^ >>yyj w 

3 



where -B(^) is a polynomial of ^. This yields that 

+C{6 + l){e-q<j)yf + Ks) + + eir)-i. 

Thus, multiplying (I4.8ip by and using the Cauchy inequality and Lemmas 4.1-4.4, 
we get 




+C(5 + 7) / / '^\G,\'d^dy + C{6 + ^){e'^Uyf + K,) ^^"^^^ 



+Cel\\i<l>yyXyy)f + C5eHl + eh)-l. 



Similarly, 




+C^(5 + 7) y J ^|Gi|^« + C(5 + 7)5^ y y ^IGJ^^ec^?/ (4.83) 

|a'|=l 
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Finally, we estimate the estimate on the highest order derivatives, that is, / e^'ijjiy(f)yydy 
and eljj^-m^d^dy with |a| = 2 in fl4.80p . To do so, it is sufficient to study the 



estimate for e / / ^-^^d^dy {\a\ = 2) because of (I4A0|) - (|4A3|) . For this, from (13161) we 
have 

vfr - ujy + Uv = e-^^vQ{f, /) = v[LmG + £^Q(G, G)]. 
Applying 5" {\a\ = 2) to the above equation gives 

= -d'^vfr + d-u.fy - J2 F""'^5"7r - 9"-"'ni9"7,] (4.84) 

|a'|=l 

+ [9"(^;LmG) - vLmO'^G] + eh'^ivQiG, G)]. 
Multiplying fITOD by eg = + yields 

v\d^f\\ 3 9"G 

= ^^{- d'^^fr + d'^uify - E - 5°-"'wi9"7j (4.85) 

\a'\=l 

+ [a°(i;LMG) - vLmO'^G] + £^9"[t;Q(G, G)]| + £t;LM9"G ■ ^ + (■ ■ ■),. 
We can compute that 



e I I \d"vfr^\d^dy 




<./l.".l/(|M.K.i|ai)^!M±^,e,, 

<C(5 + 7)e^||5°(0,^,C)f + ^^^ J J ^^^^^d^dy 



and 



\a'\=l 



eJ2 I I Id'^-^'vd-'fr^ld^dy 



\a'\=l 




|a'|=l 

a 3 /• /'^;|9°Gi2 



3 
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Similar estimates can be got to the terms ^d'^uify^^ and e d"^ °'' uid"' fy^jj^ . 

Also, we have 

(9"(wLmG) - vLmO'^G = (a"t;)LMG + 2vQ{d"M, G) 

+ J2 {2w<5(5"""'M,a"'G) + 9"-"'?;[LMa"'G + 2g(a"'M,G)]|, 

\a'\=l 

and 

e^d''[vQ{G,G)] = e^d''v)Q{G,G) + ehvQ{d''G,G) 
+ {wQ(5"""'G,9"'G) + 2(a"-"'t;)g(a"'G,G)]}. 

|a'|=l 

We only compute one of the above terms as follows, the other terms can be calculated 
similarly. 

' J J M. "^"^ 

- 16 J J ^ _ 

+Cfe5 (1 + T5r)-i + C(i + if si II (4>,. ''I',. Qf-- 

Now we estimate the term e I J vIimc^G ■ — d^dy in (I4.N,^I) . Firstly, note that 
Pi{d"M) does not contain tire term d"{v^u^O) for \a\ — 2. Tims, we have 



< -£2 




M_ J J M 

|a'|=l 

(4.86) 

Also we can get 




/ / '-W^'VWy 



^M, ^ - 16 

|a'|=l 

(4.87) 
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where the small constant rjo is defined in Lemma 4.2. The combination of (14.861) and 
(I4.87P gives the estimation oi e / vLmO^G ■ ^ d^dy. 

Thus integrating (I4.85P and recalling all the above estimates imply 




g \> . d^dy)r + -£5 




|q!'|=1 |q|=2 

+C{6 + ^)e 



\a'\=l 




(4. 



+C(5 + 7)e2 




Gi\''d^dy. 



Klel),^ ,2 



By (I4.10I) - (I4.12I) . we can choose suitable constants Cj > 1, i = 1, 2, 3, 4 so that 



|a'|=l 



>c[\\{<P,^^J,Of + e\\{K^^JyXy)\\ 
\dV\' 



d^dy 

, , 2M. ^ ^ 

l"l=2_ 

' + / / ^d^dy + e J2 

\a'\=l 





-d^dy 



l«l=2 




d^dy -C6e{l + e^T) 



1 , _ 3 
2 



Let 



\a'\=l \a\=2 



|a'|=l"^ * \a\=2 



Klel),^ ,2 




Klel),«„^,2 



(4.89) 
Gi\'' d^dy 

(4.90) 



Then by the estimates fH72|) . fOOD . (02]) . flMD . fITOD . we obtain 

+ < Cfe5 (1 + eir)-i + C56-^ (1 + £^r)-i II {%, W^j,) f . (4.91) 



5 The proof of Theorem 3.1 

Choose a large constant and set 

Eq = E2 + C^E,, Ke = K2 + C^K,. 



(5.1) 
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By combining fl4.50p and fl4.9ip . we have 

< Cofe5(i + ehy^E2 + Co6e^l + eh)-^ \\ ^y, Wy)f + Cofe^ (i + eh) 

Then Gronwall inequahty imphes that 

Eeir) < C(E6(0) + 6){1 + e^r)^, / K^iy, s)ds < C{Ee{0) + 6){1 + eh)^ 

Jo 

Now multiplying (14.911) by (1 + eh) gives 

[{l + eh)E5]r <{l+e-2T)E5r + e-^E5 

< C6e"Hl + eh)-^ + C6e'^\\{%,^y,Wy)f + e^E^ 

< C5£^l + £h)-^ +CKfi. 



Integrating (15 ■4p with respect to r and using (15.31) yield that 

E,iT)<CiEe{0) + 5){l+eh)-"K 

Thus, we have 

II VI/, W) < C\\ ($, VI/, W) II II ($„ VI/,, Wy) II < CeIeI < CiE.iO) + 5), 

and 

\\{^,^|;,C)r + e\\{<Py,^Py,Cy)r + I j^didy 

\a' r^\2 /• /• I / 1 2 




M, ^ J J m ^ ^ 

|a'|=l'' " |a|=2 

<C(E6(0) + 5)(l + £^r)-^ 

<C(i5;6(0)+5). 

And this closes the a priori estimate (14. ip . 

Now it remains to prove the decay rate of ( I3.18p . By ( 15. 5p . we have 

£^ll(0,^,C)||iso < Ce^||(0,^,C)||||(0„^„C.)l| < C7(E6(0) + 5), 

and 

1/^12 i 

' rferfz/)' + ||(^„^,)|ll(£^6(0) + 5)^ 




<C(E6(0) + (5), 
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Finally, 

/•|/(y,r,0-M[,,,,,](y,r,Or, 
sup / z-r ■ 



y 



^ ^ / m;^ + T y m:^^ (5.6) 



y 



<Ce\mij,()\\L + Cesnp J ^rf^ 

<CiEe{0) + 6)eK 
which gives flS.lSp . And this completes the proof of Theorem 3.1. 
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